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The paper reviews a concept of induced spin-flop domain inside vortices in 
an antiferromagnetic superconductor. Such phenomenon may occur when an 
external magnetic field is strong enough to flip over magnetic moments in the 
core of the vortex from their ground state configuration. The formation of 
the domain structure inside vortices modifies the surface energy barrier of the 
superconductor. During this process the entrance of the flux is stopped and a 
newly created state exhibits perfect shielding. Such behavior should be visible 
as a plateau on the dependence of flux density as a function of the external 
magnetic field. The end of the plateau determines the critical field, which 
has been called the second critical field for flux penetration. Moreover, it is 
predicted and described how this phenomenon modifies flux creep in layered 
superconductors. The various scenarios of changing the creep regime from 
thermal to quantum and vice versa at constant temperature are discussed. 
PACS numbers: i4.6O.-w, 74. 25. Ha. 



1. INTRODUCTION 

The discoveries of ternary Rare Earth (RE) Chevrel Phases REMogSs 
and RERI14B4 El compounds with regular distribution of localized magnetic 
moments of RE atoms have proved conclusively the coexistence of various 
types of magnetism with superconductivity. Intensive experimental and the- 
oretical works have shown that 4f electrons of RE atoms responsible for 
magnetism and 4d electrons of molybdenum chalcogenide or rhodium boride 
clusters responsible for superconductivity are spatially separated and there- 
fore their interaction is weak. In many of these systems superconductivity 
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coexists rather easily with antiferromagnetic order, where usually the Neel 
temperature T/v is lower than the critical temperature for superconductivity 
T c . On the other hand, ferromagnetism and superconductivity cannot coex- 
ist in bulk samples with realistic parameters. Quite often the ferromagnetic 
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order is transformed into a spiral or domain-like structure 
type and strength of magnetic anisotropy in the system 
decades the problem of the interaction between magnetism and superconduc- 
tivity has been overshadowed by high temperature superconductivity (HTS) 
found in copper oxides. However, the recent discovery of the presence of mag- 
netic order in Ru-based superconductors BHoE] has triggered a new series of 
experiments and inspired a return to the so-called coexistence phenomenon.il 
Most recently, the interplay between m agn etism and superconductivity was 
studied in d-electron UGe2i and ZrZnJH, where itinerant ferromagnetism 
may coexist with superconductivity, and in heavy fermion UPd2Al3,tLH where 
magnetic excitons are present in the superconducting phase. Among classic 
magnetic superconductors, the Chevrel phases have been studied most in- 
tensively. These compounds are mainly polycrystalline materials. However, 
some specific features can be measured only on single crystals. One such ef- 
fect is a two-step flux penetration process, predicted in Ref.(12,13|) and later 
observed solely in the antiferromagnetic superconductor (bet) ErRli4B4.ii 
This very interesting phenomenon was recently rediscovered in DyMogSs,^ 
although good quality single crystals of the classic antiferromagnetic super- 
conductors have been a long time available and measured. 

The DyMogSs compound with T c = 1.6 K exhibits transition from the 
paramagnetic to the antiferromagnetic state at Tjv = 0.4 K. Its crystal 
structure can be described as interconnected MogSs units and Dy ions. One 
such unit is a slightly deformed cube where S atoms sit at the corners and 
Mo atoms are situated at the cube-faces. The MogSs units are arranged in a 
simple rhombohedral lattice and Dy ions are located in the center of the unit 
cell. The magnetic moments of Dy ions form a simple structure consisting of 
(100) planes with moments of 8.7 [±b alternately parallel and antiparallel to 
the [111] rhombohedral axis. Neutron experiments performed on DyMogSs 
in an applied magnetic field at T = 0.2 K have revealed in the intensity 
spectrum a number of peaks characteristic for ferromagnetic order. ta These 
peaks begin to develop at Hq = 200 Oe, much below the superconducting 
upper critical field H C 2. Thus, in DyMoeSs a kind of ferromagnetic order 
coexists with superconductivity in the same manner as antiferromagnetism. 
For a field applied parallel to the [111] direction (magnetic easy-axis direc- 
tion), the ferromagnetic order is a spin- flop typeJiHI 

This feature is easy to understand. Consider the well known phase 
diagram of a two-sublattice antiferromagnet. An infinitesimal magnetic field 
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Fig. 1. On the left panel: phase diagram of two-sublattices uniaxial antifer- 
romagnet. The external field is directed along the anisotropy axis. On the 
right panel: the magnetic structure of an Abrikosov vortex. Gray area cor- 
responds to the vortex core where the spin-fjo p transition originates. This 
model has also been used by other authors 



applied perpendicular to the easy axis makes the ground antiferromagnetic 
configuration unstable against the phase transformation to the canted phase. 
On the contrary, if the magnetic field is applied parallel to the easy axis the 
antiferromagnetic (AF) phase is stable up to the thermodynamic critical 
field Ht as is seen on Fig.|l] (left panel). When the field is further increased, 
a spin-flop (SF) phase develops in the system. Let us assume that in the 
antiferromagnetic superconductor the lower critical field fulfills the relation 
H c i < \Ht and that the external field Hq is applied parallel to the easy 
axis. When H c \ < Hq < \Ht the Abrikosov vortices appear entirely in the 
AF phase. When Hq is increased beyond \Ht the phase transition to the SF 
phase originates in the core p _because near H c \ the field intensity in the core 
is approximately twice H c \ The spatial distribution of the field across the 
vortex is a function decreasing from the center as is shown in Fig. p] (right 
panel). Thus, the magnetic field intensity outside the core is less then Ht 
and, therefore, the rest of the vortex remains in the AF phase. The radius of 
a spin-flop domain grows as the field is increased. The formation of domains 
inside the vortices should be accompanied by the modification of the surface 
energy barrier. This process leads to a state of the superconductor in 
which flux entrance is temporarily prohibited, flux density in the sample is 
constant as the applied field increases (see Fig. |2|). In order to kill this state 
the external field should be increased above certain second critical field for 
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Fig. 2. Magnetic induction for DyMogSs single crystal in the virgin state 
measured as a function of an applied field for three temperatures below 
T/v = 0.4 K. The field direction is oriented parallel to the magnetic easy axis 
of the crystal. Each B{H) curve exhibits characteristic plateau indicating 
that a number of vortices is kept constant when the external field is increased. 
The calculated values of the second critical field for flux penetration H en 2(B) 
are in very good agreement with the measured onesJ^H 



flux penetration H en 2(B). Then, the vortices penetrating the sample will 
have the spin flop domains created along the cores. 

The above considerations apply to the classical superconducting Chevrel 
phases as well as to the high T c superconductors. The present work is in- 
spired by the above described discovery and the hope that the same behavior 
could possibly be observed in some of the layered superconducting structures. 
Indeed, the situation seems to be very similar in layered HTS. Here magnetic 
order is produced by the regular lattice of RE ions occupying isolating layers 
electrically isolated from the superconducting Cu-0 planes. Therefore spin 
interaction between the local magnetic moments and the conduction elec- 
trons is to weak to inhibit superconductivity. The typical example of the 
layered system is ErBa2Cri307. This compound has t etragonal unit cell with 
small orthorhombic distortion in the a-b plane r 2 H 23 l The Er ions form two 
sublattice antiferromagnetic structure of magnetic moments lying parallel 
and antiparallel to the_6 direction in the ab planeS Recently discovered RE 
nickel boride-carbidescj may serve as an another example of layered mag- 
netic superconductors. The layered structure of RE nickel boride-carbides 
is reminiscent of that of HTS and consists of RE-carbon layers separated 
by Ni 2 B 2 sheets. 0031 For example in ErN^E^C the antiferromagnetic 
structure is associated with magnetic moments of Er +3 ions, which order be- 
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Fig. 3. Schematic drawing of a piece of the layered superconductor. The 
shaded areas (n,n+l) represent superconducting layers. The bold arrows 
represent magnetic moments of RE ions lying in the isolating layers. The 
axes of the reference frame are shown. 



low 6 K in a transversely polarized planar sinusoidal structure propagating 
along a or b axis with Er moments parallel to the a or b axis respectively. 



2. LONDON THEORY 

In the following we consider the structure shown on Fig. || that we 
believe simulates a real structure of many antiferromagnetic layered super- 
conductors. A good candidate to show the above described phenomenon 
should possess the isolating layers with the magnetic moments of RE ions 
running parallel and antiparallel to the direction (easy axis) lying in the ab 
plane. 

We start description of our problem in terms of the Lawrence-Doniach 
energy functional. In this approach a layered superconductor is described 
by the superconducting planes with the interlayer distance d, as shown on 
Fig. ^. The antiferromagnetic subsystem consisting of RE ions is confined 
to the isolating layers. The magnetic moments are running parallel and 
antiparallel to the x-axis (easy axis). The Lawrence-Doniach functional is 
obtained from the standard Ginzburg-Landau energy by discretization of the 
kinetic energy in the z-direction. 

F S = [d'rdY: 
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The quantity e, m, denote the charge of the free electron and the mass of the 
current carrier in the ab plane, whereas M denotes the mass of the current 
carrier in the z-direction; A = (-4(2); Az) denotes the vector potential. The 
antiferromagnetic two sublattices subsystem with single ion anisotropy is 
described with the following energy density functional 

2 2 

f M = J2{jM ln -M 2n + K^iM^f - E (VM^) 2 } (2) 

n i=l i=l j=x,y,z 

where M n = Mi n + M.2 n is the sum of the magnetization vectors of the sub- 
lattices in the n-th insulating layer, Mf n is the component of the magnetiza- 
tion sublattice vector along the anisotropy axis in the n-th layer, J denotes 
the exchange constant between two sublattices, K is the single ion anisotropy 
constant, v^bT is the magnetic stiffness length, and Mq = |Mi n | = |M2 n |- 
The components of the total magnetization vector M have the following 
form in both sublattices :Mi y = Mosinflj, Mj x = Mq cos 6i, where 9i (canted 
spin angle) is the angle between the magnetization in the sublattice and the 
external magnetic field directed along the x-axis. The AF (6*i = 0, 62 = vr) 
and SF phases (6\ = —62 = 0) are in thermodynamic equilibrium in an 
applied field equal to the thermodynamic critical field 

H T = M [K(J-K)] 1 / 2 . (3) 

The canted spin angle of the SF phase is then expressed as 

a KM ° (a \ 

cos6 = ——. (4) 

Finally we add the magnetic field energy to obtain the free energy of the 
entire system 

F = F S + j{f M + ^(b - Mf}dV. (5) 

According to experiments the antiferromagnetic order is very weak affected 
by the presence of superconductivity, then it is reasonable to neglect the 
effect of superconductivity on the exchange interaction in F. Instead we 
introduce electromagnetic coupling between the magnetic subsystem and 
superconducting current j s . This means that both order parameters ^> n and 
M are coupled through the vector potential A 



b = rotA 



(6) 
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Fig. 4. Single Josephson vortex lying in the ab plane along the x-axis. The 
SF domain induced along the phase core is shown in the gray area. 



Mojs = rot(b — M), (7) 

where b is the vector of the microscopic magnetic field. 
The functional (§) 

can be treated in the London approximation by as- 
suming a constant modulus ^f n within the planes and allowing only for 
the phase degree of freedom. The equilibrium conditions of the system 
are the result of minimization of the Gibbs free energy functional G = 

F ' — J b(b — M)c£V. The London equations resulting from @ are following: 



d_ 

'dz' 



d 

b x + X 2 C — rot,(b - M) - A^rot y (b - M) = tpo6(y)6(z) 
d 



,d_ 

'dx 



rot z (b-M) = 







b z + A^,— rot y (b - M) - \i b — rot x (b - M) = 



dy 



(8) 



^ab = \ c y/M/m, rj = dy/M/m. The London model of continuous super- 
conductors may be used at length scales larger than the coherence length 
£, i.e. the core dimension. The structure of a vortex lying in the ab plane 
in a layered superconductor with Josephson coupling between adjacent lay- 
ers resembles the Abrikosov's one except that the order parameter does not 
vanish anywhere.Ell Instead there exists a region, rj along the plane and d 
perpendicular to it, where the Josephson current j z is of the order of the 
critical current. In this region, which plays the role of the vortex core, the 
London model fails. Away from the core the streamlines of the shielding 
supercurrents, which also represents contours of constant magnetic field, are 
elliptical except for the zigzags due to the intervening insulating layers (see 
Fig. |). 

The equations (||) should be supplemented by the appropriate set of 
differential equations describing the spatial distribution of magnetization. 
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Simple conjecture can make the calculations less complex and at the same 
time does not oversimplify the problem. We assume ™ that the magnetic 
moment is constant across the magnetic domain 



|M| = { *f % P<Pm (9) 

if p > p m 



where p m is the dimensionless radius of the magnetic domain in the coor- 
dinate system of the elliptical cylinder (x = x, y = A c pcos tp, z = A a &psin</j). 
Then the solution of Eq. @ in the cylindrical reference frame for a single 
Josephson vortex is given by the modified Bessel functions Kq and Iq 

hF = C±K (p) + C 2 Io (p) for pj < p< p m 

bAF = C 3 K Q (p) for p> p m (10) 

(pj denotes the dimensionless phase coherence length ) with the following 
boundary conditions: 

bsF (Pm) = Po#T + M = B t 

bAF (Pm) = PoH T (11) 

These conditions, together with the flux quantization condition, are used to 
calculate the arbitrary constants in (|l0|). 



BtPuJi (p, 

C 



Kx{pm) <P0 

PO-nrPr? 



_ K (p m ) 27r\ c \ ab 

p m Ki (p m ) I (p m ) - I (p m ) + p m K (p m ) h (p 

Kl (pm) <P0 



B T [pmKx (p m ) - 1] + 

C 2 = - 



PqHtp 



Kq (p m ) 27T\ c \ ab 



Iq (Pm) 

Kq (pr, 



PmKi (p m ) Iq (p m ) - J (p m ) + p m K (p m ) I\ (Pm) 

Cs = (12) 
Finally we write free energy of the isolated vortex 



A c A a 6 



f da {[bsp (r) — M] x rotbsF (r)} 



2po J o x 

+ I da {b AF (r) x rotb AF (r)} (13) 

where r = ( - — , — I is the position of the vortex line, o\ denotes the surface 



,A a 6 A c 

of the phase core, and a 2 the surface of the SF domain respectively. The 



Antiferromagnetic superconductors 



integrals in (|T^) performed as line integrals along the contours of the cross 
sections of the appropriate surfaces give s\ - the line tension of the vortex. 
The minimum of £\ with respect to p m determines 



p 2 = 50Q ^ 
m 87rA c A a fei?T 



2.1. Equilibrium energy of the vortex lattice 

The London equations can be rewritten for the lattice of vortices in the 
following way: 

B + rotrotB = -^-^5(r-r m ) (15) 
A c A ab m 

where r m specify the positions of the phase cores of the vortices. The solution 



of Eq. (15) is then a superposition 

B ( r ) = ^2 h m(r - r m ) 

m 

of the solutions b m (r — r m ) of isolated vortices at the points r m . The free 
energy of the system can thus be written as 

F = / da(B x rotB) (16) 

2^o Jo 

The above symbolic surface integral is taken over the surfaces of the phase 
cores and the surfaces of the SF domains. The energy of the Meissner state 
is chosen as zero of the energy scale. Again, when the surface integrals are 
replaced by contour ones we get line energy of the system. This, in turn, 
multiplied by vortex density n gives / the free energy density of the system. 
After some transformations one can derive the following formula 



f = ne 1 +n ( p H T (\nP)- 1 J2Ko(r m ) ; = J ^jf^ (17) 

here the sum is over all vortices excluding the one in the origin, and r m 
denotes the distance of a vortex from the origin. The lattice sum may now 
be replaced by integral in the yz-plane over a smoothed vortex density, 
excluding the area n _1 associated with the single flux line in the origin. The 
free energy density then reduces to 
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here c = |ci| denotes the length of the basal vector of the nonequilateral 
triangular unit cell, and 2|c2 1 = cV 1 + tan 2 a ( a is the angle between both 
vectors), tana = ^/j^.il To determine the equilibrium state it is nec- 
essary to minimize the Gibbs free energy density with respect to magnetic 
induction. The result yields an implicit equation for the constitutive relation 
between magnetic induction and thermodynamic magnetic field. 
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3. TWO STEP FLUX PENETRATION 

Consider a semi-infinite specimen in the half space y > 0, the vortex 
and the external magnetic field Hq running parallel to the surface in the 
x direction. The presence of a surface of the superconductor leads to a 
distortion of the field and current of any vortex located within a distance of 
the order of penetration depth from the surface. To fulfill the requirement 
that the currents cannot flow across the surface of the superconductor we 
need to introduce an image vortex with the vorticity opposite to the real 
one. Both vortices, direct and image, interact as real ones except that the 
interaction is attractive. In the low flux density regime Clem E3 has shown 
that there exist two regions: a vortex- free region of the width yff near 
the surface of the sample, and a constant flux density region for y > yff. 
Within the vortex-free area one can introduce the locally averaged magnetic 
field Bm which is a linear superposition of the Meissner screening field, the 
averaged direct vortices flux density exponentially decreasing towards the 
surface from its interior value B at y = yff, and averaged image vortices 
flux density. In our problem the x component of this superposition can be 
approximated by 

B M = B cosh { ^ ff x ~ (20) 

The boundary condition Bm(0) = HoHo determines the thickness of the 
vortex-free region 

-i ( PqHq" 



Vff = Kb cosh y— — j (21) 

We assume that the test vortex line is lying within vortex free region at a 
point r = (- — ,0j, and its image at r = (— -^-,0) outside the super- 

\A a & J ^ \ X a b ) 

conductor. Now the local field of the test vortex can be understood as a 
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superposition of the following fields 



Bsf 
Baf 

Vff 



b S F (r) - b A F (2r) + xB M (r ff - r) 
t>AF (r) - b A F (2r) + xB M {rff - r) 



(22) 



where rtf = , and x denotes the unit vector in the x direction. 

V Kb J 

Having determined the local magnetic field we can write the Gibbs free 
energy of the test vortex line as 



G 



+ 



+ 



/ da {[B SF (r) - 2/i H - M] x rotB SF (r)} 

2/i 



af (T") - 2/i H ] x rotB A F (r)} 
i da {xBm {rff — r) x rotBAF (?")} 

J cro 



(23) 



After some transformations one can obtain the Gibbs free energy per 
unit length Q 

G = Gi + G[ + G M (24) 



where 



Qi 
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Gm 
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and 
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db S F{p) 



dp 



p=Pj 



p=Pj 



dp 
db S F{p) 



db A F{p) 



dp 



p=p„ 



p=p„ 



•2pr. 



dp 
dbAF{p) 



p=Pr, 



dp 



(26) 



p=pn 



Gi describes the interaction of the test vortex with its image, G[ is a correc- 
tion term introduced by ClemJH and Gm describes the interaction energy of 
the test vortex with the modified Meissner field. To find the conditions of 
the vortex entrance and exit, one has to solve a force balance equation for 
the test vortex, at the surface of the sample, and at the edge of the flux-filled 
area, respectively. A calculation using Gi and Gm alone gives non vanishing 
force on the test vortex at r = rff. However, the force should be zero there, 
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because Gm is supposed to account for all the image vortices. To avoid 
double counting the image vortex one can subtract from the self-energy a 
contribution of the excess image fixed at r = —Vft. One can easily check 
that G[ is negligible at the surface of the sample and has no influence on the 
conditions of the flux entrance. When the flux starts to enter the sample, 
Hq = H en 2{B), 

\ v,-i ( ^H en2 {B)\ 

Vff = Ven = Aab COsh J (27) 

and the energy barrier is moved toward the surface within p m . Thus, one 
can derive from the force balance equation 



Di db AF (p) 



2D 2 dp 



P=Pm 



Ssinhf^) (28) 



The left hand side of the above equation gives H en2 (0) = HtI3{2\\i (3) 
This field may be thought as the second critical field for flux penetration cal- 
culated in the single vortex approximation c3. Combining Eqs. (^7|) and ([~~ 
one can finally obtain 



In the opposite case, when the flux exits the sample, the surface energy 
barrier tends to the edge of the flux-filled zone. Similar considerations as 
the above show that 

p H ex2 (B) ~ B (30) 
The measure of the height of the energy barrier against flux entrance is 

AH en (B) = \H en2 (B)-H eq (B)\, 

and against flux exit 

AH ex (B) = \H eq (B)-H ex2 (B)\, 

where H eq is given by Eq. ([l9|). 

Let us make a short summary of the calculations and visualize the results 
on schematic magnetization curve shown in the Fig. |5|. When the external 
field is not strong enough to create the SF domains inside vortices, than 
the magnetization process of the sample being entirely in the AF phase is as 
follows. The vortices without magnetic structure start to enter the specimen 
at H en \. When the field is increased up to the value H p i, which is of the 
order of Ht, the SF domains are created. Now, the screening current must 
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Fig. 5. Schematic drawing of the magnetization process. H en \ denotes the 
first penetration field for vortices without magnetic structure. H p i is the 
applied field which originates SF transitions inside vortices, and B p i is the 
corresponding flux density. H en 2 is the entrance field for the vortices pos- 
sessing magnetic structure. The lower dashed line shows M(H) dependence. 
The region of plateau on B{H) corresponds to the region of the second 
negative slope on M(H). 



redistribute its flow in order to keep constant the flux carried by the vortex. 
This feature is easily seen from Eqs. (|l0|)-(|l2|). The redistribution of the 
screening current changes the surface energy barrier preventing vortices from 
entering the sample as expressed in ( p8|) . It means that the density of vortices 
n is kept constant. Consequently the averaged flux density in the sample 
B = nipo remains constant when the external field is increased. In Fig. || 
this feature is visible as a plateau on the B(H) curve, or alternatively as a 
second negative slope on the M(H) curve. The vortices start to penetrate 
the sample when the external field reaches the right edge of the plateau. We 
call this value, given by J29|) , second critical field for flux penetration H en 2- 
To find the thermodynamic critical field Ht, and then to calculate 
H e n2{B) the following argumentation is used. At low fields, in the vicin- 
ity of thfi lower critical field H c \ , the intensity of the field in the vortex core 
is 2H C \ When the external field is increased the field intensity in the vortex 
core increases because of the superposition of the fields of the surrounding 
vortices. The field intensity in the core must reach Ht in order to originate 
a transition to the SF phase. Thus, taking into account only the nearest 
neighbors we can write for the nonunilateral triangular lattice 



H T = 2H cl + z ^° 




(31) 
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here c corresponds to the value B p i of the flux density for which the penetra- 
tion process stops, see Fig. |5[ From the relation B p i = 2ip ^ X a b/{c 2 \/3X c ) 
one can compute c, which in turn may be inserted back into Eq. (]3l|). It is 
easy to estimate the saturation magnetization Mq taking into account the 
volume of the elementary cell. Then, equations @ and (|j) can be used to 
calculate M in the SF-phase domain 

2KM 2 

M = 2M cos 6 = - (32) 



4. INTRINSIC PINNING 

The first quantitative approach toward intrinsic pinning in layered su- 
perconductors was based on the observation that the superconducting order 
parameter should have a periodic spatial variation across the layers. For the 
present considerations, however, the method of critical nucleus developed 
in E3 is much more convenient. The activated nucleus consists of a kink- 
antikink excitation, that is, a vortex line segment is thrown to the adjacent 
layer, thereby creating two pancake vortices of opposite vorticity, as shown 
in the Fig. [6|. The activation energy can be regarded as the energy barrier for 




intrinsic pinning. Depending on the magnitude of the driving current density 
the process may continue as the single vortex activation or the activation of 
the vortex bundle. 
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4.1. Thermal activation of vortices 

First, consider the activation of a segment of a vortex to the neighboring 
interlayer spacing (see Fig. ^). The energy associated with this process can 
be written as: 

U = 5E + V k ,-k(R) ~ U ~ jo)<PodR. (33) 

The subscript "a" or "b" of U, 5E and jo indicates that this quantity is 
calculated for a or b direction in the plane. 5E is the amount of condensation 
and magnetic domain energy that is lost at two points of the layer threaded 
by the kinks separated by a distance R. Vk,-k{R) is the kink-antikink 
interaction energy. The term proportional to the driving current j is due to 
the Lorentz force. The term proportional to jo is the energy associated with 
the distortion of the line due to the formation of nucleus. This term can be 
estimated from the simple considerations 

jwo d ~\J dydzC{y, z ) (-^ 

where the Fourier transform of the compression modulus is given by ^ 

B 2 

C(ky, k z 



a»o(i + Kt k z + 

By taking dydz ~ (po/B,u z ~ d, J| ~ k z ~ fc y (A c /A a h) the integral can be 
estimated as follows 

= (34) 

4A afe 

However, in the a direction, according to ([l4|) we have an additional contri- 
bution from the magnetic domain dydz ~ 5i^o/8-Bt- Thus we get 

5dB T , . 

joa =job+ . 2 • ( 35 ) 
U * X ab 

As the current density j drops below jo a single-vortex line can no longer 
be activated due to the confinement energy provided by the vortex lattice. 
The current jo plays very important role in all above calculations. For the 
following values d ~ 10 -9 m, £ c ~ 3 x 10 -10 m, £ a & ~ 3 x 10~ 9 m, rj ~ 10 -8 m 
and their typical temperature dependence we can estimate 



jgt„g^„io-3 B / T - 

jGL 9^0 \ T c — T 

where joL = 4eo / '{foCab^V^) is the Ginzburg-Landau critical current den- 
sity and eo = (po/(16ir 2 fJ-o^ib)- Since the depairing current is of the order of 
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10 13 [T C /(T C - T)] 3/2 then the current j ob is of the order of W W B [T C /(T C - T)]. 
Although there are no precise measurements of SF transition in the AF HTS, 
we assume that hoHt ~ 40mT. The typical value of 5.5fiB per RE atom 
per unit cell gives M ~ 0.37T. It is possible now to estimate the change of 
j'o due to the creation of SF domain along the vortex : 

J — ~ 1 + 0.625— ~3 

ha B 

providing that Ht/H c \ — 1 << 1. The energy 5E is calculated from Eq. (JsJ) 
with the right-hand sides representing the vortex cores £3 



, R 

\ x \ > 2">y = 







R 

±-,y = 0,0 < z < 



-d 



, R 

\A < —,y = 0,z 



The solution is then substituted to the free energy functional (|5|). Taking 
the limit R — » oo we exclude the energy of the kink-antikink interaction. 
Because the calculations are involved we write down only the results. 



r; 



SE b = 2de ln-^ 



(36) 



where e a = ge In 
was calculated in^ 



where 



SE a 
V k ,-k{R) 



de a In -y- 



(37) 



The energy of kink-antikink interaction 

(38) 



d 2 e /R 



(Xc/R) - Hrj/U) for 
2(\ C /R) 3 exp(-R/\ C ) for 



rj « R « A c 
R » A r 



If we introduce the quantity I ajb = 2{j - j 0a)b )/{j GL Sy/S), then Eq. (| 
be rewritten in the following way 



U„ 



f r R 
2de Un-?- + Ib— 

L Sab Sab 

r R 
dUa ln-^ + 2e I a — 

Sab Sab 



4A a / \\ c 
de / R_ 
2A a b \A C 



can 



(39) 



The critical size of the nucleus R c is given as a minimum of Eq. ( |39| ) with 
respect to R. In the approximation rj « R « A c corresponding to the 



current regime £, c d/\ 2 ab « I a b « £c/d we get 
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r>2 tl d 

AL ca,b Sao 



ab 



K = d{e a \n 




For the opposite case R » A c and £, c d/\ 2 ab » I a ,b 



Rca,b — A c In I "2 



dtic 

a,b^ab, 



U c b = 2de jln(^ ) • ; - '-)\ (J.:i) 




Kbh , 








Aab^a , / 






v - f aA a6 _ 



^ = 

When the driving current drops below jo the critical nucleus is 3D object. 
In our case it is a parallelepiped of the height R along the bundle and of the 
section S across it. The activation energy is a sum of the volume energy due 
to the Lorentz force and the surface energy. 

(BS\ 
+ j 0a ,bdR^fBS^ (42) 
fo ) 

The second term is the loss of condensation energy (and magnetic domain 
energy in the case of a direction) on both surfaces perpendicular to the 
bundle multiplied by the number of vortices threading these surfaces. The 
third term is the elastic energy released in the surface parallel to the shifted 
vortex joafidRfo multiplied by the number of shifted vortices ^/BS/ipo (one 
vortex per plane). The critical nucleus is then S c = (<£o/ B)(jo a ,b/j) 2 j Rc = 
5E a}b /(jdtpo), and the activation energy is 

K b = SE a J J ^-) 2 . (43) 



4.2. Thermal creep 



The resistive mechanism in the mixed state is determined by the acti- 
vation process leading to magnetic flux motion (creep). This motion induces 
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electric field which can be observed on the current-voltage characteristic. We 
consider the motion of activated kinks along the layers of the length L along 
the magnetic field direction. Assume that each double-kink can reach the 
boundary of the sample before the new one is created. The mean electric 
field associated with this motion is given by 



where P is the activation probability per unit volume and unit time. For 
thermal activation this probability is given by P ~ exp (—Uc/ksT) . There 
is however a crossover temperature To below which quantum tunneling of 
vortices is dominating. The probability for quantum tunneling is finite 
even for T = 0. The Neel temperature for layered antiferromagnetic super- 
conductors varies from hundreds of mK (0.6K for Erl^CusOy) to several 
Kelvin (6.8K for ErN^E^C) and therefore both mechanisms of activation 
are present in these compounds. The preexponential factors and To can- 
not be calculated in the framework of thermodynamic considerations alone. 
Fortunately, it was shown in that the activation probability of macro- 
scopic quantum excitations is proportional to j 3 . Thus we can assume that 
P = «oj 3 exp (— f7 c /fesT). Now we can calculate the current-voltage charac- 
teristics for the current density j « jo 



This almost linear dependence of E on j indicates that the resistive mecha- 
nism of bundle activation follows Ohm law. We can also calculate the rate 
of flux creep due to the thermal activation of vortices. To do this consider 
hollow cylindrical sample of a radius r and the wall thickness I « r placed 
in the magnetic field B ex > B c \ applied parallel to the cylinder axis. The 
sample has the trapped field Bi n inside the hole and corresponding trapped 
flux $ = {B in ~ B 

ex)^ r ■ According to the Faraday's law electric field due 
to the change of the trapped flux is equal to (no/2)lr(dj /dt). Combining 
this result with Eq. (El) we have finally 



E = BPLdS c 



(44) 




(45) 




(46) 



This equation can be solved analytically only in the case of the weak currents. 
Consequently for excitations in the form of bundle of vortices the above 
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equation is written as 

Ojexp<|-^(^) } + ^ = (17) 




where f2 = ^o«o Jo /(Mo7) an d 7 = rl/(Ld) is the factor determined by 
the geometry of the sample. The solution of Eq (|47| ) is given in terms 
of exponential integrals and for the case of jo a ,b/j — 1 << 1 it can be 
approximated as: 

2 

| \n(l+uj a , b t) (48) 



j(0) , $(0) 1 fc B T /j(0) 



j(t) $(t) 28E a>b \j 0a , b 



4ip a 5E ab ( j 0a ,b\ J 5E ab ( j 0ab 

Wa,b = Ti-rr exp 



where 

_ 4^ 

M07 Vi(0);^\ fc 6 T Vi(0). 

This result is in agreement with the experiments on HTS.H For << t << 
1/uj the change of trapped flux is linear in time and for t >> 1/uj logarithmic. 
In antiferromagnetic superconductors, however, we see additional change of 
characteristic frequency as the magnetic field changes its direction in the ab 
plane. 



4.3. Quantum creep 

Below To, the quantum activation probability is essentially indepen- 
dent of temperature P ~ exp (—S/h) and is interpreted as arising from the 
quantum tunneling of vortices through intrinsic pinning potential.E^'Ed In 
the folowing we show a considerable change of tunneling rate and crossover 
temperature due to the SF phase transition around the vortex core. Now, 
consider the vortex line as a straight string-like object of an effective mass m 
per unit length trapped into a metastable state in an intrinsic pinning poten- 
tial V{u) and exposed to continuous deformation u(x, t) in the z direction.0 
The magnetic field is applied in x direction (a direction on Fig. In 
the semi-classical approximation the quantum decay rate is calculated as a 
saddle-point solution (bounce) of the Euclidean action S for the string 



S 




1 fdu\ 2 £l (du\ 2 T _, , 



, 2 

u(x, t) — u(x, T ) 



2vr 



r — r 



(49) 
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Here (3 = (ksT)^ 1 , 77 is the viscosity coefficient and r denotes imaginary 
time. The pinning potential V(u) consists of intrinsic periodic part and the 
Lorentz potential: 



V{u) 



fojpd 
2vr 



■ cos 



2ttu 



ipoju. 



(50) 



For large current, this potential can be expanded around the inflection point 
to give 

' u \ " f U 
w) \w , 



V (u) = V Q 



(51) 



where Vn 



*wfl w 3 and w 



'o - 3 

of the barrier because V(0) = 
called Caldeira-Leggett action 



¥ ( 



Jo -J 
2j 



may be thought as the width 



w) = 0. The last term in (49) is the so- 
which describes ohmic damping produced 
by the coupling to the heat-bath of harmonic oscillators. The line tension 
£\ is different for vortices in two different orientations in the ab plane. As 
previously discussed the vortices lying parallel to b direction and those laying 
in the a direction but created in the magnetic field fulfilling relation H c \ < 
H < \Ht S have the line tension equal to 



60 d ' 



(52) 



For those vortices lying in the a direction but possessing spin flop domain, 
we write the following expression^] 



ip H T 9 , tpo 

£o In — ot^. 

2 128 vrrfS T 



(53) 



In the semiclassical approximation the decay rate is given by the value of the 
action on a classical trajectory obtained from the Euler-Lagrange equations 
of the motion 

d 2 U d 2 U i , , 77 [ h P du 7T / /\ 

- m d^- £l d^ + v ^ + npl dT ^ cot np( T - T )=° (54) 



The trajectory uq(x) for static solution of Eq. (|54|) gives the activation energy 
in the thermal regime T > Tq. Below this crossover temperature a new kind 
of trajectory, periodic in imaginary time, develops. Therefore, u(x, r) can 
be expanded in the Fourier series with Matsubara frequencies 



U(X, T) = ^2u n (x) COS (0J n T) 
n=0 



27T77, 



(55) 
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Substituting this expansion into (|5^ ) and linearizing potential around the 
static solution uq(x) one obtains 

d 2 u n 



dx 2 



+ V"(u )u n = - (nu) n - mu^j u n . (56) 



The above equation has three discrete solutions,c3 the unstable one corre- 
sponding to the tunnelling process determines the crossover temperature 



Airm 



20Trip j m f j - j\2 



drj 2 V 2j 



(57) 



The above calculations apply to both kind of vortices. The only difference is 
their effective mass and viscosity coefficient. It is possible to express these 
parameters as the function of condensation energy accumulated in the vortex 
cores. For the stationary flux flow the viscous force is equal to Lorentz 
force. The electric field generated by the moving vortex is E = B^, so we 
get E = ^^-j = pj = PNjf^j where pn is the normal phase resistivity in 
the ab plane and is the upper critical field parallel to the layers. Finally, 

<PQ H c2 VoKH c y/2 4\/3k 2 

rj= = =ei , (58) 

PN PN np N In K 

where H r = £ ' k2 1 v ^ is calculated from the constitutive relation e-i = H r \^n. 

L TTLpQ In K PH 

The effective mass of the vortex can be deduced from the work of Suhl.eil 



3_ eni 



He derived the core contribution to the inertial maSS Tn C ore — g Mt e — \ ; 

where m e denotes the mass of the electron and is the Fermi energy, and 
the electromagnetic contribution coming from the energy of the electric field 
induced by the moving flux. Simple estimation shows that this contribution 
in layered superconductors is 10 -4 of the core contribution. Therefore, 



2 _9>^gnePo_ 
(p^1T 2 €F (In k) 

It is possible now to relate the crossover temperature in Eq. ([57]) to the line 
tension of the vortex 

T = ae^ 1 . (60) 
The coefficient a depends on the material constants and current intensity. 



4.4. Crossover from quantum to thermal creep 



As was already mentioned there are two types of vortex lines in the 
system. The first ones, without magnetic domain, occur when the field 



T 
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thermal creep 
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quantum 
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77- 
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Fig. 7. Schematic diagrams showing possible ways ( marked by arrows ) 
of changing quantum creep behavior in the system to thermal one, and vice 
versa. Shaded areas on the diagrams correspond to quantum creep regime eJ. 



is applied in the a or 6 direction, but its intensity does not exceed \Ht- 
Eq. gives their line tension and the related crossover temperature T b ~ 
er 1 . The other type, possessing magnetic domain, occur when the field is 
applied in the a direction and its intensity exceeds \Ht- Eq. (|53| ) gives their 
line tension and the crossover temperature as T^ a ~ e" 1 . It is easy to see that 
e a > Eb and therefore T$b > To a . The above calculations lead to the following 
conclusion. It is possible to switch the creep regime at constant temperature. 
To do this, one needs to change the field intensity or simply change the field 
direction in the ab plane. The diagrams in Fig. |7| show possible scenarios of 
crossover from quantum to thermal regimes. Let us discuss just two of them. 
The first prescription is following. Fix the temperature To somewhere in the 
range T f, a > Tq > To a . Then align the external field in the a direction and 
increase its intensity to the point marked " 1" on the left diagram in Fig. [?]. 
The system is in the quantum creep regime now. Then increase the external 
field beyond \Ht- The system jumps to the point "2" of the left diagram 
and finds itself in the thermal creep regime. Doing the same operations 
in the reverse order one enforces the system to crossover from thermal to 
quantum creep regime. The other scenario is the following. Apply magnetic 
field along a axis and increase its intensity above \Ht keeping temperature 
constant in the interval Tq\, > Tq > To a . Then move the direction of the 
external field from a to b axis. The system goes now from point "2" of the 
left diagram (thermal creep) to the point "3" of the right diagram (quantum 
creep) . 
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